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Abstract 

A recent variant of the inflationary paradigm is that the "primor- 
dial" curvature perturbations come from quantum fluctuations of a 
scalar field, subdominant and effectively massless during inflation, 
called the "curvaton" , instead of the fluctuations of the inflaton field. 
We consider the situation where the primordial curvature perturba- 
tions generated by the quantum fluctuations of an inflaton and of a 
curvaton field are of the same order of magnitude. We compute the 
total curvature perturbation and its spectrum in this case and we 
discuss the observational consequences. 

1 Introduction 

The first observations of the WMAP satellite have confirmed the 
basic predictions of the inflation scenarios [2] (see also ) , namely the 
existence of super-Hubble primordial fluctuations, mostly adiabatic 
and characterized by a quasi-scale-invariant spectrum. Because direct 
information on the primordial universe is scarce, it is essential, when 
trying to interpret the cosmological data, to determine to which extent 
a precise measurement of the primordial spectrum gives us information 
on the inflationary scenario itself. 

In this perspective, whereas in the standard inflationary paradigm, 
the primordial perturbations are generated by quantum fluctuations 
of the inflaton, it has been realized recently that the observed pri- 
mordial fluctuations could be accounted for by quantum fluctuations 
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of a scalar field other than the inflaton, dubbed the curvaton field 
^ El (and see also jTj and [S]). The curvaton is a scalar field a 
whose effective mass must be much smaller than the Hubble param- 
eter during inflation so that it acquires quantum fluctuations with a 
quasi scale-invariant spectrum. After inflation, the primordial curva- 
ture perturbations are generated by the conversion of the isocurvature 
perturbations due to the fluctuations of the subdominant curvaton 
field, into adiabatic perturbations. 

In the curvaton scenario, inflation is still necessary since this is the 
source of the curvaton quantum fluctuations but the separation of the 
roles played by the inflaton that drives the evolution of the universe 
and by the curvaton which generates the perturbations, liberates the 
inflaton from some observational constraints. This alternative possi- 
bility is mixed news for the inflationary scenarios. It is good news in 
the sense that it allows to construct more easily explicit inflationary 
models with a firmer footing from a particle physicist's point of view 
[SJ. But it is also bad news in the sense that it shows a new limita- 
tion to extract information on the nature of the inflaton field from the 
measurement of the primordial fluctuations. 

In the curvaton scenario, the curvature fluctuations generated by 
the inflaton are usually assumed to be negligible, and the primordial 
fluctuations are thus essentially due to the curvaton. However, as 
was pointed out in , one can envisage cases where the fluctuations 
generated by both the inflaton and a curvaton-like field are relevant. 
This is the situation that we explore in the present work. 

It can occur, for instance, when is of the order of the Planck 
mass mp during inflation, still assuming that the curvaton energy den- 
sity is much smaller than that of the inflaton (otherwise we would be 
in the context of double inflation ^J), and implies that the curvaton 
starts to oscillate when its energy density already contributes signifi- 
cantly to the total energy density. This is in contrast with the usual 
curvaton scenario where the curvaton starts oscillating long before it 
dominates and can thus be treated as a pressureless fluid during the 
first radiation dominated phase. 

Curvaton and inflaton generated perturbations can also be of the 
same order of magnitude for small a^, i.e. cr* <C mp, but with a slow- 
roll inflation parameter e sufficiently small during inflation to produce 
big inflaton perturbations. 

This paper is organized as follows. In the next section we write the 
background evolution equations and we derive their initial conditions 
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explaining the context of our present analysis. In Sec. Ill we discuss 
the perturbation equations and their initial conditions. We can an- 
alytically derive the curvature perturbation after curvaton decay, i.e. 
at the onset of the standard big bang cosmology, in two limits. In 
the pure curvaton limit, for cr^, ^ mp, we recover the curvature per- 
turbation already derived in previous works. In the limit where the 
curvaton acts as a secondary inflaton, for cr* > mp, we recover the 
perturbation of the double inflation model. In the intermediate case 
we resort to a numerical analysis to derive the mixed inflaton and cur- 
vaton perturbation, which is the main result of our work. In Sec. IV 
we derive the spectrum of the curvature perturbation, its spectral in- 
dex, and discuss the observational consequences of the model. Finally 
we conclude in the last section. 



2 Homogeneous equations of motion 

In the curvaton scenario, the inflationary phase is followed by a first 
radiation dominated era, where the dominant species are the decay 
products of the inflaton. The curvaton is a scalar field, a, that is 
assumed to be subdominant during the whole inflationary epoch and 
at the beginning of the radiation dominated post-inflationary era. Its 
effective mass is assumed to be much smaller than the Hubble pa- 
rameter during inflation. Here, for simplicity, we take the quadratic 
potential 

V{a) = \m'a\ (1) 

In the post-inflationary era, two components coexist: the radiation 
fluid and the curvaton. The evolution of the homogeneous background 
is thus governed by the Friedmann equation 

H' = '-f {Pr + + ^mV^) , (2) 
the continuity equation for the radiation energy density pr, 

Pr + AHpr = 0, (3) 
and the equation of motion for the curvaton, 

u + 3H& + Tr?a = 0. (4) 
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When H ^ m the cosmological friction is so strong that the cur- 
vaton is essentially frozen. This is the case during inflation since 
the curvaton must be effectively massless to acquire the usual scale- 
invariant quantum fluctuations. After inflation, H keeps decreasing 
until it reaches the value ~ m, at which point the curvaton starts 
oscillating. 

In the standard curvaton scenario, a is supposed to start oscillating 
while it is still strongly subdominant with respect to radiation. Using 
(PI, this implies 

~ » ^al (5) 

TJlp 

where we have introduced the reduced Planck mass mp = (SvrG)"^/'^. 
This relation means that the initial amplitude of the curvaton must 
satisfy cr# <^ mp. Conversely, assuming an initial value a* ~ mp 
during inflation implies that the epochs of oscillation and domination 
roughly coincide for the curvaton. 

During the radiation dominated era, when a is subdominant, H = 
(2t)~^, and the curvaton equation of motion @ reads 

3 

a + —a + m^a = {). (6) 
One recognizes a Bessel equation, whose non-decaying solution is given 

by 

{mtYl^' "2V4r(3/4)' 

where Jjy stands for a Bessel function of rank v. In the limit mt <C 1, 
one gets (T(j) ~ a^,. The expression (T(j) given above embodies the 
initial conditions for the curvaton field. 

Note that the solution Q differs from the slow-roll solution that 
would apply during inflation when the Hubble parameter is slowly 
varying. Here, the acceleration is not small with respect to the other 
terms in the equation of motion In particular, in the limit mt <C 1, 
one finds 



-H'^(j) =^ -^mV(j). (8) 
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3 Analysis of the perturbations 
3.1 Equations of motion 

The equations governing the evolution of the scalar-type perturbations 
are well-known and can be found in review papers such as or jl^]- 
We adopt here the longitudinal gauge in which the perturbed metric 
reads 

ds2 ^ ^ 2<^)dt^ + a^{t){l - 2^)6ijdx'dxK (9) 

The perturbation of the curvaton field is denoted 5a and the per- 
turbations of the radiation fluid are described by the perturbation of 
the energy density, 5pr-, or the density contrast 5r = 5pr/ Pr-, and the 
peculiar velocity Vr- 

Taking into account the relation ^' = imposed by the Einstein 
equations in absence of anisotropic stress, the perturbed equations 
of motion are the energy and momentum constraints in Einstein's 
equations, respectively 

- 3if(i?$ - — $ = AtiG [dSa + w?aba - a^^ ^ pr^r] (10) 



and 



6 + = 47rG {aba - ^PrV^ , (11) 
the continuity and Euler equations for the radiation fluid, respectively 



^;,-4$ = (12) 
and 

Vr - HVr + <!> + -6r = 0, (13) 

and finally the Klein-Gordon equation for the curvaton, 

Sa + 3H6a + + = ^o"^ - 2mV$. (14) 

Since the cosmological perturbations of observational interest to- 
day were on super-Hubble scales in the early universe, we will be 
interested only in the long wavelength limit k/{aH) — > of the above 
equations, which reduce to the following system: 

- 2,H{H^ + 6) = 47rG (aSa + n?a5a - + pr5r) , (15) 
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5r-A^ = 0, (16) 

and 

5a + 3H6(J + m'^6a = - 2n?a^, (17) 
where we have gotten rid of the perturbed velocity Vr- 

3.2 Initial conditions 

Our initial conditions are defined early in the first radiation domi- 
nated era when the curvaton is at rest and subdominant. In this 
limit the Bardeen potential <I> is constant. In the standard curvaton 
case, it is further assumed to be zero, but here it cannot be neglected 
since it represents the contribution of the inflaton to the "primordial" 
curvature perturbation. Our initial conditions are thus given by 

$ = $^, $ = 0, 6(T = 5(T^, 5ct = 0, (18) 

where <I>* stands for the curvature perturbation generated during infla- 
tion. Using the energy constraint H15|) . these initial conditions imply 

6r = -6— (5fT*. (19) 

Pr Pr 

Moreover, since the curvaton is initially subdominant, the second term 
on the right hand side can be dropped and the initial condition for 6r 
reduces to 

5» = -2$,. (20) 

The evolution of 6a is governed by the equation of motion Eq. ()17() . 
which can be approximated in our limit by 

3 ■ 

6a H 6a + m^Ja = -2m^$*o-. (21) 

One recognizes the same Bessel equation as for the background but 
now with a source term on the right hand side, which is proportional 
to the background solution for a given in Eq. The solution for 
the perturbation is then found to be given by 



(22) 
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Note also that using Eq. we can rewrite Eq. (|^^ as 



6a, 



6a^, 



{iy- 



(23) 



This yields (5cT(j) = (5ct^, + which justifies the initial conditions 

for the curvaton perturbation given in (|18|) . 

As the energy density perturbation for the curvaton is given ex- 
plicitly by 

6pa = rm^aSa + aSa - fT^^>, (24) 
we can compute, using (|^ . the initial density contrast 

. (25) 



= 2— - - — 
cr* 2 



3.3 Evolution of the perturbations 

The above analytical expressions are valid only during the early ra- 
diation dominated era when the curvaton contribution to the back- 
ground energy density is negligible. This contribution however keeps 
increasing with time until it eventually becomes significant, and a 
more general analysis is required. 

The purpose of this section is to compute the curvature perturba- 
tion after the curvaton decay, i.e. to establish the initial conditions at 
the start of the second radiation dominated era, which must be identi- 
fied with the usual radiation era of the standard cosmological scenario. 
Since we assume the decay to be instantaneous (it has been shown re- 
cently that this is a very good approximation [TT), it is enough for 
our purpose to determine the curvature perturbation when the cur- 
vaton is both dominating and oscillating, which will be denoted <!>/. 
The corresponding curvature perturbation in the subsequent radiation 
phase is then simply obtained by applying the usual transfer coeffi- 
cient between a matter dominated phase and a radiation dominated 
phase. 

In order to determine the final curvature perturbation in our sce- 
nario, i.e. the primordial perturbation for the standard cosmological 
model, we have solved numerically the system of equations that con- 
sists of ©-(EI for the background and of (fT3)) - H17() for the perturba- 
tions. For the homogeneous system, the only parameter that can be 
varied is the initial amplitude cr* of the curvaton. As for the perturba- 
tions, the initial conditions depend on two parameters, the curvature 
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perturbation and the curvaton perturbation (5(T*. Since this is a 
hnear problem, the final result can be written, quite generally as 



= a^^ + /36a^, (26) 

where the two coefficients a and /3, to be determined, depend only 
on the homogeneous quantities. This means that one can analyze 
independently the cases <I>^= = and 5a^ = 0, the general result being 
given by the sum of these particular cases. 

Let us start with the case da^, = 0. We will now show that this 
case corresponds to a purely adiabatic initial condition. Let us first 
recall that the so-called isocurvature perturbation between two matter 
components (here radiation and curvaton) can be defined as (see e.g. 

EI) 

5 = 3(Cr-Ca), (27) 

where for a given matter species X, is the curvature perturbation 
on uniform energy density hyper surf aces, namely 

PX 3{l + wx) 

with wx = Px/px- For the curvaton, the pressure is = {d"^ — 
m^(T^)/2 and the energy density is = (<3"^ + m?a'^)/2 so that the 
entropy perturbation is given by 

- %.-%S.. (29) 



1 + Wr 1 + Wa 4 (T^ 

Substituting the initial conditions (PU]) and ()25() , one gets for the initial 
entropy perturbation 

S« = -2^^. (30) 

Therefore, Jc* = corresponds to S*^*^ = 0, i.e. to a purely "adiabatic" 
initial perturbation. 

The perturbation then remains adiabatic throughout the evolution 
from the radiation dominated era to the oscillating curvaton dom- 
inated era. Note that, in contrast with the case of two uncoupled 
adiabatic perfect fluids where the individual (x are separately con- 
served, the entropy perturbation does not in general remain constant 
on super-Hubble scales because a non zero intrinsic entropy of the 
scalar field implies Co- 7^ 0. However, if the initial entropy vanishes it 



8 



will remain so subsequently. Since we have a pure adiabatic pertur- 
bation, the total C is conserved on super-Hubble scales. Consequently 
the curvature perturbation is simply governed by the global equation 
of state of the matter content, or in other words the total ( is con- 
served in time. During a cosmological era with w = const, the Bardeen 
potential $ is constant and is related to C via the relation (see e.g. 

m) 

C = -4^<f. (31) 
^ 3(l + '«;) ^ ^ 

Using this relation both in the initial radiation dominated era with 
w = 1/3 and in the oscillating curvaton dominated era with w = 0, 
the conservation of C implies that the coefficient a in Eq. ()26() is given 
by 

a = —. (32) 
10 ^ ^ 

We have checked numerically that this result indeed holds for all types 
of initial conditions (not only Jc* = 0). 

To determine the other coefficient, /?, it is now sufficient to consider 
the initial configurations with = 0. We distinguish below several 
cases: first the limiting case where a behaves as a standard curvaton, 
i.e. oscillates long before domination (see upper part of Fig. 1); then 
the other limiting case where a starts oscillating only after it com- 
pletely dominates and thus generates a second inflationary phase (see 
upper part of Fig. 2); and flnally the general intermediate case, as 
illustrated in the upper part of Fig. 3. 

3.3.1 The standard curvaton limit 

This is the limit corresponding to <^ inp. As discussed before, 
this implies that the curvaton oscillates long before its domination. 
One thus recovers the standard curvaton scenario, in which case the 
curvaton can be identified with a pressureless fiuid during the post- 
inflation phase. 

Let us briefly recall the usual analysis in the standard curvaton 
scenario jlj. The curvaton is assimilated to a dust-like fluid with 
density contrast 5m- Since the two fluids are decoupled, Cr and Cm are 
separately conserved, whereas the total ( is given by 

^ ^ 4prCr + 3pmCm ^^^^ 
4pr + 3pm 
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Figure 1: The standard curvaton limit, for o"* = 0.3mp. Evolution of the 
curvaton, its energy density fraction fl„ = Pa/{Pa + Pr) (upper part), and 
of the curvature perturbation $ (lower part). The initial conditions for the 
perturbations are Sa* — mp and = —2. 
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Assuming no initial curvature perturbation, i.e. ^>=k = and thus = 
0, one finds that during the curvaton domination 



C = Cm = ^(^«, (34) 

where the second equahty comes from the definition of ^ given in 
Eq. (|28j) . Using the relation (|31|) between and <I> during aw = era, 
one gets 

^f = -lc = -lsS. (35) 

The final step is to relate the density contrast to the initial curvaton 
fluctuation, by using Eq. (j^^)- For = 0, this is simply = 
25a^/a<t:. Therefore, in the standard curvaton limit, one finds for the 
coefficient (3 of Eq. the value 

P = -^. (36) 

OCT* 

Although the expression for a has already been obtained by the 
general argument given earlier, it is nevertheless instructive to see 
how this value is recovered in the curvaton limit by considering a non 
vanishing <I>^,. When averaged over several oscillations, Eq (|25|) yields 

^ 5« = 2^ - . (37) 

(J* z 

During the curvaton dominated era, one finds (note the difference with 
Eq. (El) 

C = Cm = -^* + ^5^). (38) 

Substituting in = — (3/5)C, derived from pTjl . one finally gets the 
expected result 

9 2 5cr* 

We have also checked this result by evolving numerically the pertur- 
bation equations. As an illustration, we show the evolution of <I> in 
Fig. 1 for initial perturbations 6a^ = mp and = —2 ^. 



^Within perturbation theory, the overall amplitude of the perturbations can be ar- 
bitrarily rescaled and we have chosen to set Sa^ ~ mp. Moreover, the initial condition 

= —2 corresponds to the curvature perturbation produced by a chaotic quadratic model 
with 0* = 6mp and 6(j)t = mp. 
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Note finally that the above equation can be interpreted as a par- 
ticular case of the more general equation 

^/ = + (40) 

relating the final curvature perturbation, in a matter dominated era, 
to the initial curvature perturbation and entropy perturbation defined 
in a radiation dominated era. Eq. H40|) can be easily derived from 
$j = — (3/5)Cm by substituting Cm = Cr — S/3, which follows from 
(Ei, with Cr = + = -(3/2)$,. The expression ()4U() can 

be used only when the entropy perturbation is conserved on large 
scales. This is the case for the curvaton only when it is oscillating or 
when it dominates (and its intrinsic entropy thus vanishes), because 
otherwise Co- is not conserved. In the early radiation dominated era, 
the expression for the entropy perturbation is given in Eq. I|3()j) . When 
a oscillates, Pajo'^ = 1 and S is constant with the value S = — 2(5(T*/(T*. 
Substituting in (|4()j) one indeed recovers (jHOI)- 



3.3.2 The secondary inflaton limit 

Another limiting case corresponds to the situation where a starts os- 
cillating only after it completely dominates the universe. This implies 
that the curvaton a is responsible for an additional phase of inflation: 
this is why we call a the secondary inflaton in this situation. 

We now attempt to evaluate analytically the final curvature per- 
turbation. Let us consider again the entropy perturbation S. As men- 
tioned already, S is not constant because Co- is not conserved. This 
can be checked explicitly with the expression 

5 = _2^^, (41) 

already given in Eq. (|30|) . Indeed, in the early radiation dominated 
era, using the explicit solution given in ((71), one finds the behavior 
p^/c3-2 ~ (25/8)(mt)-2. 

In the subsequent phases dominated by the curvaton, i.e. the sec- 
ondary inflation and the oscillating phase, the entropy S remains con- 
stant. Although the expression (|41|) is rigorously valid only during the 
phase when the curvaton is negligible, we can evaluate it at the tran- 
sition between the (first) radiation dominated era and the curvaton 
dominated era to get an estimate of the subsequently constant value 
of S. 
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Figure 2: The secondary inflaton limit, for a^, = Snip. Evolution of o", 
its energy density fraction (upper part), and of the curvature perturbation 
(lower part). The initial conditions for the perturbations are Sa^ = m-p and 
= -2. 
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During the curvaton dominated phase, a fohows a slow-roll motion, 
characterized by 



2 2 9 

^ o ma ma , , , 



and therefore, the coefficient Pal^'^ takes the constant value 

Pa 3 crj 



^ - 7^' (43) 
0-2 4 mf, 

where we have used the fact that a is essentially frozen at its infla- 
tionary value cr* throughout the radiation dominated era. Inserting 
this ratio in Eq. ||1T|) . we get 



S = --A^ba,. (44) 
imp 

Substituting the above value in Eq. (|4()|) , we then find for the curvature 
perturbation in the strongly oscillating phase 

= _ (45) 

^ 10 10m| ^ ^ 

We have checked numerically that our analytical derivation gives a 
very good approximation. This is illustrated in Fig. 3 for two different 
values of (and with ba^ = mp). 



3.3.3 The general case 

We now consider the general case, which includes the intermediate 
situations where the curvaton starts oscillating at about the same 
time it becomes dominant. We have computed the dependence of the 
coefficient /3 on the initial value o"*, by solving numerically the full 
system of equations for different values of a^,: the intermediate case 
(T* = mp is plotted in Fig. 3 as an illustration. 

The resulting function, up to a normalization constant defined just 
now, is plotted in Fig. 4. Instead of using /?, we prefer to work with 
the coefficient that appears in the expression for the curvature per- 
turbation after the decay of the curvaton, i.e. in the second radiation 
phase of our scenario, identified with the radiation dominated era of 
the standard big bang model. The relation between the curvature 
perturbation at the end of the curvaton dominated phase and the 
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2 4 6 8 10 12 

N 

Figure 3: The intermediate case, for = mp. Evolution of cr, its energy 
density fraction (upper part), and of the curvature perturbation (lower part) 
for two different initial conditions: = and for $^ = —2 (both cases with 
5(7* = mp). 
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Figure 4: The function f{a^), which characterizes the amphtude of the con- 
tribution of (7 to the curvature perturbation. For cr* <^ mp one recognizes 
the oc 1/(7* contribution of the pure curvaton modeL For cr* ^ mp, one 
recognizes the contribution of a secondary infiaton, proportional to cr* . 
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curvaton perturbation $_r_d in the post-curvaton radiation phase is 
given by 

the transfer coefficient being simply due to the transition from a mat- 
ter dominated phase to a radiation dominated phase. It is then con- 
venient to define the function / as 

^RD = - ^Sa, . (47) 
rrip 

As we can see on Fig. 4, for small values of a^,, we are in the standard 
curvaton regime with 

4 Trip 

f{(ytf) = , pure curvaton, (48) 

9 cr* 

whereas for large values of cr=K, we recover the secondary inflaton limit 
characterized by 

/(cr*) = — secondary inflaton. (49) 
3mp 

In the intermediate range, / reaches its minimum value, / ~ 0.55 
around cr* ~ 1.2mp. 



4 Observational consequences 

4.1 Power spectrum of the primordial fluctu- 
ations 

So far, we have computed the primordial fluctuations, defined in the 
second radiation era after the decay of the curvaton, without specify- 
ing the initial curvature perturbation However, in the context of 
an inflationary era driven by a single slow-rolling scalar field, one can 
relate to the fluctuation of the inflaton. The perturbations due to 
the inflaton can be computed following the standard technique for a 
slow- rolling single field inflation model (see e.g. JS])- One finds 

54)^, (50) 



$. = -[l-(l + 3C)6 + C^]-^^ 
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where V = V{(f)) is the potential of the inflaton, V' = dV/d(j), e and rj 
are the first two slow-roll parameters, 



2 W 



(51) 



2^" 

V = mp — , (52) 

and C = —2 + ln2 + 7 ~ —0.73, 7 being the Euler constant. It is 
important to include the slow-roll corrections in the normalization of 
the curvature perturbation since these can be of the same order as the 
curvaton contribution. The subscript * means that all quantities have 
to be evaluated at Hubble crossing during inflation. Note that we are 
implicitly assuming that large scales that are observable today, i.e., 
those of order Hq^ , did not reenter the Hubble radius between the 
end of inflation and the curvaton domination. 

Substituting in our result of the previous section, Eq. (|37j), the 
primordial curvature fluctuations in the second radiation era are thus 
given by the combination 



$«z. = -[l-(l + 3C)e + C,]^^ 



/((T*)5(7*. (53) 

mp 



The fluctuations 50* and 6(t^ must be seen as independent random 
fields with the same amplitude for their power spectra and can be 
written 

H* = TT^e^, (5fT* = —^ea, (54) 
zvr zvr 

with (e^ecr) = 0. It is also worth noticing that our result, in the 
secondary inflaton limit discussed before, with /(<t*) = a^/{3mp) is 
in agreement with the expression obtained in (see also ^H]) for 
the curvature perturbation given in terms of the perturbations of the 
two scalar fields in a double inflation model. 

Let us drop the subscript *. Using (|53|) and (|5l|) . the amplitude of 
the power spectrum for <I> is 

V^ = [l + pe - 2(1 + 3C)6 + 2Cv] ,,^J ^ , (55) 



where we have introduced 



ISvr^emp 



/ = 4/- (56) 
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The contribution of the curvaton to the total curvature perturbation 
appears here on the same footing as the slow-roll corrections in the 
amplitude of the inflaton power spectrum. However, when / S> 1 the 
curvaton contribution dominates over the slow-roll corrections and the 
term —2(1 -|- 3C)e -|- 2Crj can be neglected. 

The relative importance of the curvaton over the inflaton pertur- 
bation is parametrized by 

^ ^ l-2(l + 3C)e + 2C7? ^^^^ 

This is given by the product of two parameters that depend on two a 
priori independent physical processes: e, which must be smaller than 
one, is determined by the inflationary scenario, while p > 1.36 is de- 
termined by the curvaton expectation value. The slow-roll parameter 
e is typically of order ~ 1/60 in chaotic models of inflation but it can 
be much smaller in other inflationary models. Therefore, there are 
cases where the curvaton value during inflation is small, a <C m-p, and 
yet the perturbations generated by the inflaton are of the order of or 
larger than the curvaton perturbations, i.e. R <j 1. 

When the curvaton perturbation is absent, i.e. for / = 0, we re- 
cover the pure inflationary spectrum (see e.g. jlSj'l. 

V^ = [l- 2(1 + 3C)e + 2C77] -^g^2-^2 ■ (58) 
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Conversely, when the curvaton perturbation dominates the power spec- 
trum, i.e. for R ^ \ and a <C mp, we recognize the spectrum of 
perturbations of the curvaton , 

= 8W ' 

where we have used / = (2^2/3) (mp/a) from Eqs. (|1H|) and 

The presence of mixed curvaton and inflaton perturbations also 
affects the scalar spectral index, which is defined as 

Inserting in this expression the spectrum (|55j) and assuming that the 
curvaton mass is much smaller than the inflaton mass, m <^ V" , one 
finds 

- 1 = -2e+ + [-2{7 + 12C)e^ + 2{3 + 8C)e7j-2Cf], (61) 

1 + f^e 
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where we have introduced the second-order slow-roll parameter 

V'V" 

e^m%^^. (62) 

In the square bracket we have isolated terms which are second order in 
the slow-roll parameters. These can be as important as the curvaton 
contribution if / ~ 1. 

Neglecting second order terms, in the case of pure inflation, = 0, 
one recovers the well-known result 

n^j — 1 = 2r/ — 6e, (63) 

whereas in the case of a pure curvaton, /^e » 1, one gets |1] 

n,-\ = -2e. (64) 

The power spectrum for gravitational waves is given by the usual 
expression (see e.g. jTJ]) 

2 rj2 

7'r = [l-2(l + C)e]^^. (65) 

Using (|55() . we thus find for the tensor-scalar ratio 

4 Vt 16e 
y 1 -I- /^e 

(we have defined r as in JZI)- Neglecting second order terms in the 
slow-roll parameters, for / = we recover r = 16e, the tensor-scalar 
ratio for a pure inflationary model. For a pure curvaton perturbation, 
/^e 3> 1, one finds r ~ 16//^ ~ 18((T=K/mp)^. In the usual curvaton 
context, where a* ^ mp, this implies that gravity waves are negligible 

El- 

Adding a curvaton perturbation to a pure inflationary perturbation 
always decreases the tensor-scalar ratio r. When / is large, the effect 
on the scalar index depends on the sign of 2r/ — 4e: if it is negative 
(this is the case for the so-called small field and chaotic models |18j). 
the scalar index is increased, which means that the spectrum is bluer; 
if it is positive, the effect is opposite. 

Models of inflation are often classified according to the value of 
their first two slow-roll parameters e and r\ ^Hl ■ Since for / = there is 
a one-to-one correspondence between the (e,r/)-plane and the {ns,r)- 
plane, it is customary to use the latter in order to easily confront 
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models with observations. In the mixed case, / 7^ 0, the presence 
of the curvaton introduces a degeneracy between these two planes, 
which depends on the parameter cr^,. In the next subsection we will 
use Eqs. (1^1]) and to show the effects of this degeneracy for the 
special case of A(/>^-inflation. 

Let us now discuss the consistency relation. We restrict the discus- 
sion to the large / case and we neglect terms which are second order 
in the slow-roll parameters. The tensor spectral index being given by 
the usual expression 

"^=^h^ = -'^' 

the consistency relation between the tensor-scalar ratio and the tensor 
spectral index is modified into 

-'"^ (68) 



1 - PnT/2 

For large / we find a modification of the standard consistency relation 
valid for a pure single-field inflation. This represents a distinctive 
feature of mixed curvaton and inflaton perturbations, which breaks 
the degeneracy of mixed models. Indeed, this modification depends 
solely on the parameter a^,, which determines the value of /: since the 
relation (|68j) does not depend on a particular inflaton potential, it can 
be seen as a way of measuring a*. 

For completeness, we briefly discuss the running of the scalar spec- 
tral index. We find 

-4.(,-2e)-l^f^^i^!^ + i^l(!^. (69) 



dink " ' l + /2e (l + /2e)2 

For / = we recover the pure inflationary running |17j . 

"^""•^ -246^ + IGer/ - 2^^ (70) 



din A: 



while for /^e 3> 1 we find the running for a pure curvaton spectral 
index, 

In the next subsection, we will illustrate the general formulas de- 
rived here by considering the particular case of an inflaton with a 
quartic potential. 
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4.2 Mixed perturbations with V{(f)) = 

The quartic potential V = A(/)^ is a simple inflationary potential which 
has attracted some attention lately as it lies in a region excluded by 
the WMAP data analysis |17l 1191 l2Uj . For this particular inflationary 
potential, we now illustrate how mixed curvaton and inflaton pertur- 
bations can be constrained by the data and how the observational 
predictions of a pure inflationary model change in the presence of a 
curvaton. 

The slow-roll parameters for this potential are 

2 2 4 

The number of e-foldings before the end of inflation when the pertur- 
bations at our present Hubble scale exited the Hubble radius during 
inflation is given by 



iV*=ln^ = -^, (73) 

a* OTTlr 



*^end 

where we have assumed that (/)end 0*- Using this expression we can 
write the slow-roll parameters in terms of the number of e-foldings 
AT,. 

Determining the appropriate number of e-foldings A^* is very im- 
portant, since it tells us which part of the potential corresponds to 
the present observable fluctuations. This depends on the history of 
the universe from the end of inflation until today. In particular, it is 
very sensitive on the duration of the reheating phase: this is usually 
the largest uncertainty in its estimation. However, in the case of pure 
inflation with a quartic potential, A'^* can be evaluated very accurately 
because the expansion during reheating is the same as in a radiation 
dominated era. Indeed, the energy density of an oscillating scalar field 
in a quartic potential behaves as radiation, and the duration of the 
reheating phase is thus no longer important. In this case, the number 
of e-foldings can be evaluated precisely and is 

^quartic ^ (74) 

This number assumes radiation domination from the end of infla- 
tion until equality. It does not take into account a possible interme- 
diate phase of curvaton domination between reheating and nucleosyn- 
thesis, which lasts 

AiV = -ln^ (75) 

3 Pdec 
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e-foldings, where pcB is the total energy density at curvaton and radi- 
ation equahty and pdec ~ SF^mp is the energy density of the universe 
when the curvaton decays with decay rate T. Curvaton domination 
thus decreases the number of e-foldings which becomes 

TV, -64- -3- In (76) 

t^ Pdec 

Let us compute this correction in more details. If ~ rrip the 
domination phase starts with a small period of inflation. When the 
curvaton and radiation energy densities are the same, we have for the 
total energy density 

PCD = 2p^ ~ m^al (77) 

In this case we obtain 

1 777 1 (T 

iV* ~ 64.1 - -In— - -In— , for a, ~ mp. (78) 
6 r 6 mp 

It is also interesting to consider the pure curvaton case correspond- 
ing to (T* <^ mp . In this case Eq. ((77|) does not apply, but we can still 
compute the ratio pdec/pcD- Once the curvaton starts oscillating at 
H m, we have 

2 2 f ^osc\'^ o 2 2 /'^osc\ 

p^~m^<(— — I , p,. ~ Sm^mf, ( — — I , (79) 

where Oqsc is the scale factor at the onset of the curvaton oscillating 
period. This implies 

acD / flosc ^ 3m| /al- (80) 

On the other hand, when the curvaton decays {H ~ F), it dominates 
the universe and 

f — V =^ SrVf, . (81) 

On combining Eqs. (|8fl|) and (|81|) we can compute the correction factor 
in Eq. ()76() . and we obtain finally 

1 Tfi 2 (T 

iV* ~ 64.4 - - In -- - - In — , for a* <. mp. (82) 
6 r 3 mp 

The ratio m/F depends on the specific model for the curvaton. 
It parameterizes the duration of the curvaton oscillating phase and 
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represents here the main uncertainty, although hmited by the prefac- 
tor 1/6. Also the last term can be important for very small cr*: the 
smaller (T^,/mp is, the later the curvaton starts dominating the uni- 
verse, thus shortening the duration of the curvaton domination phase. 
This increases the number of e-foldings and partially compensates the 
ln(m/r) correction. 

Now we have all the ingredients to study the predictions of mixed 
models as compared to a pure quartic inflationary model. We use 
Eqs. (|72() and (|73() to express the slow-roll parameters e and rj in terms 
of the number of e-foldings A^*. Furthermore, on using Eqs. and 
H66|) . we can express the observables Ug and r in terms of A*, and 
we can compare the predictions on the (ns,r)-plane of a pure quartic 
model and mixed models with cosmological data. 

In Fig. Elthe two dimensional likelihood contours at 68% and 95% 
confidence level of the WMAP data (combined with other data as from 
the analysis of fSOl) are shown on the (re<j, t) -plane. For a pure infla- 
tionary model, the = 64 realization is marginally excluded by the 
data. When adding the perturbation of the curvaton, the predictions 
on the (ns,r )-plane are changed and the A* = 64 realization of the 
inflationary model can be safely included into the 95% confidence level 
contour for a mixed model with < 0.5mp, and into the 68% con- 
fidence level contour for a mixed model with < O.lmp. However, 
one must also take into account the fact that the number of e-foldings 
N^t: decreases because of the phase of curvaton domination, accord- 
ing to (|76j) . As an illustration of this effect, we consider the interval 
60 < A^, < 64 on the curvaton trajectory and we plot it in Fig. [21 
Decreasing the number of e-foldings tends to exclude the model. 

As anticipated in the introduction, the extraction of information 
from the data is limited now by the introduction of a new field, the 
curvaton, and a corresponding new parameter, the curvaton expecta- 
tion value during inflation. The degeneracy of the data can be broken 
by making use of the consistency relation, Eq. (|68|1 . although we are 
still far from measuring n-r j22j . 

5 Conclusion 

In the present work, we have explored one particular facet of the cur- 
vaton mechanism: we have considered the case where primordial cur- 
vature perturbations generated by the curvaton are of the same order 
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Figure 5: Two-dimensional likelihood contours at 68% and 95% confidence 
level of the WMAP data on the (n^, r)-plane, as compared to the predictions 
of a pure inflationary model (solid A^-trajectory), a mixed model with / = 1 
(upper dotted A^-trajectory), corresponding to cr* ~ 0.5mp, and one with / = 
3 (lower dotted A^-trajectory), corresponding to a* ~ O.lmp, respectively. 
The stars denote A^ = 64 on the A^-trajectories. The interval 60 < A^ < 64 
is represented by bold dots on the A^-trajectories of the mixed models. The 
likelihood contours are from the analysis of S. Leach and A. Liddle pU] . 
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of magnitude as the traditional inflaton-generated curvature pertur- 
bations. This can happen when the vacuum expectation value for the 
curvaton is of the order of the Planck mass during inflation. This situ- 
ation can occur in some scenarios, such as in the string axion scenario 
discussed in 

Independently of the specific model of inflation, we have obtained 
the expression for the "primordial" curvature perturbation, given as a 
linear combination of the post-reheating curvature perturbation (i.e. 
generated during inflation) and of the curvaton perturbation. In terms 
of the curvature perturbation on the uniform energy density hyper- 
surface this is given as 



where H, a, and e are the Hubble parameter, the vacuum expectation 
value of the curvaton, and the first slow-roll parameter, all evaluated 
at Hubble crossing during inflation, respectively [see Eq. H53() ]. This 
expression is only valid for large /, otherwise corrections which are 
second order in the slow-roll parameters must be introduced. The 
function /(cr) depends on the initial background value for the curvaton 
and interpolates between two regimes, which have been investigated 
previously in the literature: the standard curvaton regime, where the 
curvaton starts oscillating long before it dominates, and for which 
/ = (2\/2/3)(mp/o"), and the secondary inflaton regime, where the 
"curvaton" is still frozen when it dominates and thus produces an 
additional period of inflation, and / = (l/\/2)(o"/mp). 

In the case of slow-roll inflation, we have computed the spectral 
amplitude and index associated with curvature perturbations of mixed 
origins and shown how the resulting expressions nicely interpolate 
between the pure inflaton and pure curvaton cases. As an illustration, 
we have considered the case of a quartic potential for the inflaton and 
shown how allowing for a mixing with a curvaton changes the usual 
predictions. 

The results derived here are quite generic and can be applied to 
any inflationary model. Although we have considered, for simplicity, 
only a quadratic potential for the curvaton, which is a good approx- 
imation near a local minimum of most potentials, it is in principle 
straightforward to extend our results to other potentials. As we have 
shown, a generic consequence of mixed perturbations is a modification 
of the usual consistency relation, which gives some hope to be able to 




r 1 



H 



(83) 
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distinguish observationally between the standard inflationary scenario 
and its mixed version. 
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